Transport with three-particle interaction 
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Starting from a point - like two - and three - particle interaction the kinetic equation is derived. 
While the drift term of the kinetic equation turns out to be determined by the known Skyrme mean 
field the collision integral appears in two - and three - particle parts. The cross section results 
from the same microscopic footing and is naturally density dependent due to the three - particle 
force. By this way no hybrid model for drift and cross section is needed for nuclear transport. 
Besides the mean field correlation energy the resulting equation of state has also a two - and three 
- particle correlation energy which are both calculated analytically for the ground state. These 
energies contribute to the equation of state and lead to an occurrence of a maximum at 3 times 
nuclear density in the total energy. 



I. INTRODUCTION 

The equation of state of nuclear matter is known to 
saturate at the density of no = 0.16fm~ 3 with a binding 
energy of — 16MeV. Refined two-particle calculations us- 
ing Bruckner theory M and beyond are not crossing 
the "Coester line" in that they lead to binding energies 
and/or densities above the saturation ones. Only den- 
sity dependent Skyrme parameterizations originally in- 
troduced by Skyrme |||| , three-particle interactions m 
or relativistic treatments can reproduce the correct 
binding energy and saturation density of nuclear matter. 
The Skyrme parameterizations are derived from an effec- 
tive three-particle interaction ||. This leads to a non- 
linear density dependence in the effective three - particle 
part which is responsible for saturation. The importance 
of three-particle collisions in nuclear matter transport has 
been demonstrated, e.g. in ||. The density dependence 
deviating from that arising by three-body contact inter- 
action has been introduced and compared with experi- 
ments in 1 1 1 . 

The relativistic approach on the other hand yields im- 
mediately the correct saturation with two - particle ex- 
change interactions. This has been traced down to the 
nonlinear density dependence of scalar density and con- 
sequently nonlinear mean fields |11| which leads to a 
density contribution to the binding energy of (n/no) 8 ' 3 . 
Higher order effects have been shown to lead to a 
(n/no) 3 ' 4 dependence, see 12 and citations therein. 



Therefore the physics of the saturation mechanism is 
certainly a nontrivial density dependence of the mean 
field. Though the Skyrme interaction has been over- 
whelmingly successful and the evidence of saturation 
by three-particle interaction m, it is puzzling that a 
transport theory with three-particle contact interaction 
has not been formulated. In this paper we will derive 
the corresponding kinetic equation using noncquilibrium 
Green's function and we will show that the three-body 
interaction term leads to a (n/no) 10 / 3 contribution to the 



binding energy. 

The transport theory including three - particle interac- 
tions seems to be of wider interest e.g. in [ [l3[ it has been 
found that the three-body interactions have a measurable 
influence on thermodynamic properties of fluids in equi- 
librium. The description of noncquilibrium in nuclear 
matter is mostly based on the Boltzmann (BUU) equa- 
tion which uses the Skyrme parameterization to deter- 
mine the mean field and drift side of the kinetic equation. 
The collisional integral is then calculated with a cross 
section which arises from different theoretical impact. In 
addition one usually neglects the three-particle collision 
integral. These hybrid models have worked quite success- 
fully despite their weak rigour in microscopic foundation. 

In this paper the transport equation will be derived 
for two- and three - particle contact interactions. Wc 
will see that a natural density dependent mean field ap- 
pears, in agreement with variational methods. Further 
we will obtain, on the same ground, density dependent 
cross sections. This has the advantage that we can de- 
rive a BUU equation which has the same microscopic 
impact on the drift and collisional side. Finally the 
three -particle collision integral appears naturally from 
this treatment. From balance equations we will derive 
the density dependent energy which gives the equation 
of state. 



II. THREE - PARTICLE KADANOFF AND BAYM 
EQUATION 

We consider a Hamilton system with the Hamiltonian 



i 



£ a~ta+a+a ia: ja k V 3 (i, j, k) 



(1) 



tjk 



where af , ai are creation and annihilation operators with 
cumulative variables i = (xi,U, ...). We assume the po- 
tential is contact - like as 



V 2 (i,j) = t 8(i-j) 
V 3 (iJ,k) = t 3 S(i-j)S(i-k). 



(2) 



The Heisenberg equations of motion for the annihilation 
operators read 

ihdt 1 ai = [<zi, H] 



h 2 V 2 

-^-ai + V V(l, 2)ata 1 a 2 
- Y V(l, 2, i)a^at,a 1 a 2 a 3 . 

1 23 



(3) 



From this we get the equation of motion for the causal 
Green function G(l, 2) = A < Ta^a 2 > with the time or- 
dering T and the average taken about the nonequilibrium 
density operator 



(ihd tl - 



)G(1, l')=5(l-l') 



2m 

d2V(l,2)G 2 (l,2,l',2+) 



where the 1 + denotes the space - time point 1 = (x\ , t% + 
e) with a time infinitesimally larger than 1. This equa- 
tion is enclosed in the standard manner introducing the 
self energy E. We have split the mean field parts Y^hf in- 
cluding exchange and have introduced for the rest the self 
energy on the Keldysh contour. This is equivalent to the 
weakening of initial correlations. The correlation func- 
tions are G<(1,2) =< a^ai > and G > (1,2) =< aiaj > 
respectively and their equation of motion follows from (^) 
in the form of the known Kadanoff and Baym equation 



i(Go 1 G < - G < G(7 1 ) = E fi G< - G<E A 



(5) 



where one takes advantage of the retarded and advanced 
functions G R ' A = =R#[±(*i - t 2 )](G> T G<) and under- 
stands products as integration over inner variables. 




FIG. 1. The definition of self energy. We denote 
the two-particle interaction by long dashed lines and the 
three-particle interaction by short dashed line. 



A. Mean field parts 

We first calculate the mean field parts. These parts 
come from the first order interaction diagrams of figure 
|l|. To this end we use the two- and three- particle Green 
function in lowest order approximation seen in figure 



- g: 



— — X" 



~ J d2d3V(l,2,3)G 3 (l,2,3,l',2 + ,3+) 
= (5(1 — l') — J dx 2 'SHF(xi,x 2 , t)G(x 2 ,t, xi,t) 
■ [ d2(E(l,2)G(2,l') -E > (1,2)G < (2,1')) (4) E HF (l,l') = 



FIG. 2. The lowest order expansion of Green functions. 
The statistical factor g due to degeneracy in the system is 
indicated explicitly. 



Introducing these diagrams into figure [I], one obtains 
for contact interaction (0) 



<?±1 



t n(l) 



(g±l)(g±2) 
2g 2 



t 3 n(l) 2 



5(1-1') 
(6) 



where we have used the fact that the density is n(l) = 
G < (1, l)/g. Here and in the following we write the up- 
per sign for Bosons and the lower sign for Fcrmions. The 
result (|^) is the known Skyrme mean field expression 
in nuclear matter for g ~ 4. It resembles an effective 
density dependent two-particle interaction arising from 
three-particle contact interaction. As a consistency check 
we see that for g = 2 degenerated Fermionic system, like 
spin -1/2 Fermions, no three-body term appears. Pauli- 
blocking forbids two particles to meet at the same point 
with equal quantum numbers which one would have for 
three - particle contact interaction and degeneracy of 
(7 = 2. 



B. Kinetic equation 



For the kinetic equation we need as the lowest order 
approximation the next diagrams of figure || including 



2 



one interaction line. While there are different expansion 
techniques at hand |l6| , |l3|| we will perform here a scheme 
as close as possible to the causal Green's function includ- 
ing all exchanges which are presented in figure |^. This 
will give us the advantage that we can consider all dia- 
grams which differ by exchange of outgoing lines as equal. 
To this aim we introduce the abbreviated symmetrized 
Green's function with respect to incoming lines in figure 




FIG. 3. The symmetrized three-body Green function with 
respect to incoming lines which are used in figure 

X = zr 



Since we can also have nontrivial combinations of two 
incoming lines and one free line for the three - body 
Green function we have introduced an auxiliary two - 
body Green function G| which is defined in figure |^. 
Please note that the exchange of outgoing lines does 
not lead to distinguished diagrams. The two - particle 
Green function is then given by figure ^. Introducing 
now this expansion into the definition of self energy in 
figure [l] we obtain the 4 diagrams of figure q. All other 
combinations lead either to disconnected diagrams or to 
equivalent ones by interchanging outgoing lines. Here the 
seemingly disconnected diagrams for G\ in figure | lead 
together with the corresponding diagrams from G\ to the 
vertex in front of G\ of figure g. 




X = 



o 



FIG. 4. The introduced vertices. The three - and two - 
particle interaction potentials are the same as in figure^]. 

The expansion of the three -body Green's function is 
given in figure |^. Besides the three -particle interaction 
potential we have to consider all possible single line in- 
teractions between two incoming lines. This is shortened 
by the introduction of the symmetrized Green function 
in figure |^ and the vertices of figure ^. 
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FIG. 5. The first order expansion of the two - and three 
particle Green's functions with respect to interaction. 





FIG. 6. The self energy according to figure^ with the ex- 
pansion of figure ^| 

The enclosing diagram about the auxiliary Green func- 
tion G| is given in figure |?]. Please denote that here the 
seemingly disconnected diagrams of figure || vanish. We 
have used again the fact that the interchange of outgoing 
lines does not lead to distinguished diagrams. 







FIG. 7. The enclosing diagram of the auxiliary Green func- 
tion G|. 

Inserting figure in figure || we obtain the final result 
for the self energy in the Born approximation 



S<(1,1') = (5±1) [tl 



+ (.9 ± 2) 4* 



6* 



,n(l)n(l') 



• 5*0*3 ( 



n(l) , n(V) 



xG < (i,yyG > {i'i) 
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-~(g±l)(g±2)t 2 3 G<(1,1') 3 G>(1'1) 



(7) 



Note that all interactions vanish for the fermionic case 
g = 1 since the Pauli-principle does not allow contact 
interaction in this case. For g = 2 only the two - parti- 
cle interactions survive for fcrmions as one expects from 
the Pauli principle. That the result is symmetric in the 
densities at the two time-space points needs no further 
comment. 

From the self energy we can write down the Kadanoff 
and Baym equation (^) in closed form. If we employ 
standard gradient expansion and convert the equation (||) 
into an equation for the pole part of the Green's function 
we obtain the final kinetic equation 



at m or 

2 f dp 2 dp[dp' 2 2 . 

t 



x / drcos 
't 



(2nh) 
1 



P2-P1 -P2) 



(E 1 + E 2 - E[ - E' 2 )(t - t) 



x {f[{T)&{T)h{T)h{T) - h(r)f 2 (r)f[(T)f^r)} 



+ 



2 f dp 2 dp 3 dp' 1 dp' 2 dp 3 2 , , , 

-T 3 (pi +p 2 +P3 ~Pi -P 2 -P3) 



(2-iTh) 12 

(E 1 + E 2 



x / drcos 

'to 

*{f[(T)ti(T)&(T)h(T)J 2 (T)h(T) 
h{T)h{T)h{T)J[{T)j' 2 {T)m} 



T? 1 J?' J?' s 

Ei~ E % - E 3i 



(8) 



with / = 1 — /, and the particle dispersion E = p 2 /2m + 
<7 HF . The quasiparticle distribution functions are nor- 
malized to the density as g J ^th)* /(p) = n - ^ or sa ke 
of simplicity we have suppressed the notation of centre 
of mass space dependence. Neglecting the retardation in 
the distribution functions and taking to — > 00 gives the 
standard Boltzmann two- and three- particle collision in- 
tegrals. 

The introduced two - and three- particle T-matrices 
are read off from (]?]) as first Born approximation 



To = 



9 2 
(9 ±2) 



iti + 6t 



,n(l)n(V) 



5t t 3 ( 



n{\) n(l') 



^a 2 = 5^(9 ±2)4 



divide by and for the two- and three - particle out- 
going collision partners respectively. An effective cross 
section can be defined in the usual way and read up to 
constants 



The reader is remind that we have generally to sum over 
collision partners and to average over the outgoing prod- 
ucts. While the summation is already taken into account 
explicitly in the diagrammatic approach above we had to 



da 2 
~dfl 



d<r 3 

dn 



T 2 = — 
2 16 



<)/- 



-tin 2 
4 3 



5t a t 3 n 



27 

Tg 



1 n2 2 . 1 . 

(to + -t 3 n) + -t 3 n(t + -t 3 n) 



T 2 



64 



(10) 



If one compares this with the effective two- particle cross 
section derived in jl7| one sees that it differs by the last 
two terms of ( fio|) proportional to t 3 . This is due to the 
fact that the underlying three - particle process is taken 
into account explicitly here while in [fl7| an effective two 
- particle kinetic theory has been developed. Obviously 
we have to face a partial cancellation between three - 
and two - particle collisions. This will become explicit in 
the discussion of correlation energies. We will see that 
the two - and three - particle correlation energies indeed 
cancel partially concerning the t 3 terms. 

The ratio of the effective cross sections are for the case 
of nuclear matter (g = 4) 



da 2 
da 3 



T 2 
T 2 



9^r + —n + 5 — n 

4 4 *3 



(11) 



This could serve as the measure of relative importance 
of the corresponding collision processes, but we prefer to 
discuss this later in terms of dispersed energy by the two- 
and three- particle collision integrals since this includes 
also the Pauli blocking. 

We would like to point out that we have fulfilled our 
task and have derived a kinetic equation where the drift 
as well as the collision integral follows from the same mi- 
croscopic impact. The drift represents a Skyrme mean 
field and the collision side shows a two - and three - par- 
ticle collision integral where the cross sections are calcu- 
lated from the same parameters as the mean field. 



III. BALANCE EQUATION 

By multiplying the kinetic equation with X,p and E, 
one obtains the balance for density n, momentum u and 
energy £. Since the collision integrals vanish for density 
and momentum balance we get the usual balance equa- 
tions 



dn d_ r _dp_dE f = 

[ ' dt OR J (2nh) 3 dp 3 

dtii . d f dp dE 

— 1 - ^ 3 (r— f + £oij) = 



dt dRj J (2-Kh) 3 ^ l d Pj 
where the mean field energy of the system varies as 



(12) 
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se = 



dp S£ 
{2nhf 5f(p) 
dp 



■Sf(p) 



such that from (|^) follows 



E5f(p) 



fc 2 

£ =< — > 
2m 



— ^ — ton 



(£±1XS±2). 3 
6? t3U - 



(13) 



(14) 



With the help of this quantity the balance of energy den- 
sity reads 



d£ d 



dt 8RJ (2irh) 3 dp 



dp dE d d 

E—f = -—E co „ 2 (t)-—E corT3 (t) 



(15) 



with the two-particle correlation energy JT^,|T^| 



-^COIT2 (^) 



~iy* 

to 



dp 1 dp 2 dp' 1 dp' 2 , , 

■T 2 6{p 1 +p 2 -p 1 -p 2 j 



(2nhf 



x sm 



(Ei + E2 — E'i — E' 2 



mf' 2 (r)h(r)h(r) 
(16) 

and the complete analogous expression for the three par- 
ticle energy. 

Please note that in the balance equations for the den- 
sity and momentum no correlated density or correlated 
flux appears. This is due to the Born approximation. 
For more nontrivial approximations as e.g. the ladder 
summation these correlated observables appear po| , pT| . 



Then the nuclear binding at tiq — 0.16 fm 3 with an 
energy of Eb = — 15.68MeV is reproduced by the set 

t = — (2e f - 5E B ) = -1026.67MeVfm 3 

15n 

16 « 
h = — (e/ - 5E B ) = 14625Mevfm 6 
5n 

which leads to a compressibility of 

K = 9n2 ^J^ = jj^n 2 , - ^e/ = 377MeV. (18) 



2. Parameterization with two - and three - particle 
correlation energy 

Now we consider the two - particle correlation energy. 
From (@) we obtain the total energy 



£ + E corn + E 



(19) 



We want to calculate the long time limit which 
represents the equilibrium value. From the identity 

oc 

/ sin xtdt = — a principle value integration replaces the 


sin term in (|16| 



Eco„ 2 (00) 



-9 



dp\dpidp' x dp' 2 
(2n7if 

x5(pi +p 2 - pi -p^fifzhh- 



To 



Ex 



J7> TP' Tpf 

£J2 ~ &\ ~ &2 



(20) 



For ground state Fermi distributions this expression 
can be integrated analytically jl9| and we find from j20| ) 
the known Galitskii result for the two - particle ground 
state correlation energy p2| 



A. Fit to nuclear matter ground state 

As one can see, it is not enough to use the mean field 
parameterization to fit the equation of state as done in 
most approaches so far. Since the collision integral in- 
duces a two - particle and three - particle correlation 
energy we have to take this into account and have to 
refit the parameter to, £3 to the binding properties of nu- 
clear matter. To illustrate this fact let us evaluate the fit 
without two - and three - particle correlation energy and 
therefore without collision integrals. 



1. Hartree-Fock parameterization 

Taking into account only Hartree-Fock mean field cor- 
relations ( |l4| ) the ground state energy for nuclear matter 
reads 



E, 



E £ 3 3 1 2 

— 7 = — = — €f H — tnn H fan 

A n 5 / 8 16 



(17) 



. . 2 log 2 - 11 f PfmT 2 . 2 
= 4e / ^ 1 , 9.3 ) 



35 
5.7910" 5 
MeVfm 2 ' 



,4/3 



(21) 



As pointed out in p3fl we had to subtract here an infinite 
value, i.e. the term proportional to /1/2 in (p(i|). This 
can be understood as renormalization of the contact po- 
tential and is formally hidden in E COIT (0) when time in- 
tegrating equation (|1J). For finite range potentials we 
have an intrinsic cut-off due to range of interaction and 
such problems do not occur. 

The correlational two - particle energy ((2^) is always 
negative for fermionic degeneracies 2 < g < 8 and for 
bosonic degeneracy g < 4. Since the leading density goes 
with n 10 / 3 it will dominate over the positive kinetic part 
which goes like n 2 / 3 . We find that for densities around 
3no the total energy has a maximum and starts to de- 
crease continuously for higher densities. The maximal 
energy is changed towards higher values if we now in- 
clude three - particle correlations since the latter remain 



5 



positive and have the same leading density term as the 
two - particle correlational energy. 
The three - particle part reads 



£corr 3 (OO) = -g 



dp 1 dp 2 dp 3 dp' 1 dp' 2 dp' 3 2 

-L<i 



{2irk) 15 3 

x 6( P1 +P2+P3- p[ -p' 2 - P^f'xfLf'Jihh 
V 

v 

Ei + E2 + E3 — E[ — E' 2 — E% 



(22) 



The analytic result for the contact potential has not been 
given in the literature to our knowledge and reads 



4e/ 



9013 



pjmT 3 ^ 

2-9-25-77-13M7r 4 ?i 6 ' 
2.3710" n io/3^2 
MeVfm 2 3 ' 



(23) 



This three - particle correlation energy remains positive 
and has the same leading density behaviour of n 10 / 3 as 
the two - particle correlational energy. Since the pre- 
factor is smaller than the one for the two - particle corre- 
lational energy we obtain a maximum at 3 times nuclear 
density beyond which the two - particle part dominates 
and the total energy diverges negatively which would 
mean a collapse of the system. This clearly marks the 
limit of the Born approximation. 

Taking now the correlational energy into account via 
( |l9| ) instead of only the Hartree Fock energy ( |l7| ) we ob- 
tain a fit to the nuclear binding of 



t = -745.71MeVfm 3 
U = 8272.8.MeVfm 6 

which leads to a compressibility of 

K = 351MeV 



(24) 



(25) 



which is somewhat lower than the mean field compress- 
ibility of (|l8|). The comparison of the two equation of 
states with and without two - and three-particle corre- 
lations can be seen in figure |[ The inclusion of two 
- particle correlation energy leads to a maximum at 3 
times nuclear density above which the system collapses. 
The complete result including three particle correlational 
energy leads to a higher reachable maximum. 
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FIG. 8. The equation of state due to mean field compared 
with the equation of state including two - and three - particle 
collisions in Born approximation. 



B. Importance of three particle collisions 

We have seen that the three - particle contact inter- 
action induces in a natural way a density dependence of 
the two - particle cross section. Additionally we have ob- 
tained an explicit three - particle collision integral on the 
same microscopic footing. We want now to answer the 
question how important the three - particle collisions are. 
Therefore we use as a measure the ratio of the two- and 
three - particle ground state correlation energies, j2l| ) 
and (p3[), since this gives the measure of how much en- 
ergy is maximally dispersed by the corresponding inte- 
grals. We obtain 



Eo 



2 5 - 5 ■ 11 • 13 
9013 



(11 -2 log 2) 



3 r t 
n 2 t{ 4 tm 



u 3 i (, 3 

This ratio is decreasing until it reaches the density 
18t 



(26) 



5t 3 



0.32fnT 



where the ratio has the minimum 
E 2 



1.4. 



For higher densities the constant value 
E 2 



E't 



18.3 



(27) 



(28) 



(29) 



is approached. 

I other words this means that the importance of three 
particle collisions increase with increasing density up to 
twice the nuclear density where the correlational energies 



G 



are nearly equal. For higher densities we have 18 times 
larger two -particle correlational energies than three - 
particle ones. At nuclear saturation density this ratio 
is 



E2 

E~3 



= 19.2 



(30) 



indicating that the three -particle collisions become im- 
portant between nuclear density and twice the nuclear 
density while it can be neglected in the other cases. 



IV. CONCLUSIONS 

For a microscopic two - and three - particle contact 
interaction consisting of two parameters the Kadanoff 
and Baym equation of motion is derived. From this a 
Boltzmann kinetic equation is obtained with drift which 
turns out to be the known Skyrme Hartree Fock expres- 
sion while the collision side consists of two - and three 

- particle collision integrals. The two - particle collision 
integral contains an explicit density - dependent cross 
section arising from the three - particle contact interac- 
tion. By this way both the drift side as well as collision 
side is derived from the same microscopic footing and 
no hybrid assumptions about separate density dependent 
mean field and cross section is needed. 

The correlational energy for the three- and the two 

- particle part are calculated analytically. Due to the 
density dependent two - particle collision integral both 
correlational energies have the same leading power of 
(n/n ) 10 ^ 3 4 in density which shows that both contribu- 
tions are of the same importance if three - particle inter- 
action has to be considered. This is clearly motivated by 
the saturation point where non - relativistic two - particle 
approaches fail to overcome the "Coester line" . 

We find that there is a maximum in the energy at 3 
times nuclear density if two - and three particle corre- 
lational energies are included. Beyond this density the 
Born approximation fails at least in that the system col- 
lapses towards diverging negative energy 

While the two - particle collision cross section has a 
natural density dependence due to three - particle contact 
interaction it turns out that the explicit three particle 
collision integral can be neglected as long as one is below 
nuclear density. Around twice nuclear density the three 
particle collision integral has the same importance as the 
two - particle one since it disperses the same amount of 
energy. For higher densities the three - particle collision 
integral is again negligible. 
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the manuscript. 
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